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Abstract 

Evolution of structure functions and their moments at low and moderate is 
studied in the chiral field theory. Evolution equations based on perturbation ex- 
pansion in the coupling constant of the eff'ective theory are derived and solved for 
the moments. The kernels of evolution arising from different processes have been 
calculated with contributions from direct and cross channels, the interference terms 
being non-negligible in the kinematic regions under consideration. This is shown to 
lead to flavor- dependence of the kernels which manifests in observable effects. The 
invalidity of the probabilistic approach to the evolution process is also pointed out. 
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1 Introduction 



While perturbative QCD has been reasonably successful in describing the scale 
dependence of the structure functions at asymptotically large Q"^ [1], phenomena 
observed experimentally at intermediate seem to demand consideration of non- 
perturbative effects[2], the observed departure of the Gottfried sum from the parton 
model value providing a notable example [3] . Effective theories based on qq'K vertex, 
which are expected to incorporate the relevant nonperturbative effects, give simple, 
and at least qualitative, explanation of the observed phenomena [4- 6]. 

In such theories, the QCD evolution of the effective parton densities due to quark- 
gluon interaction is replaced essentially by evolution due to quark pion interaction. 
The evolution equations proposed in these theories involve evolution kernels or split- 
ting functions which can be computed from the virtual photon cross sections for pion 
emission processes. While in perturbative QCD, the evolution kernels result directly 
from the Weiszacker- Williams procedure and admit simple probabilistic interpreta- 
tion, we point out in this paper that the correct evolution kernels in pion-emission 
theories are not consistent with such interpretation. 

Although we illustrate the point with the specific example of the chiral field 
theory [7] used in [4], the conclusions are absolutely general. 

The effective quark densities are defined in the Altarelli -Parisi formalism through 
the use of leading log approximation. In effective theories, the lack of asymptotoic 
freedom requires a different procedure for defining the effective quark densities. We 
discuss this in the next section. In section 3, the evolution equations for the effective 
parton densities thus defined are obtained. In section 4, we also show the evolution 
equations that would follow from the probabilistic approach in a kinematic region 
where such approach may be supposed to be valid and point out the difficulties that 
would result from such a procedure. In section 5, we discuss the evolution of the 
moments of the structure function that follow from the evolution equations derived 
in section 3. 

2 Effective quark densities 

The standard theoretical approach to the Q^-dependence of the structure functions 
begins with the formula 

Ff{xM^)lx [ ^g.o(y)<7^**(x/y,g^) (1) 

where g^o is the naive parton model density of the i-th parton in the target nucleon 
and cr'^*'^' is the ■j*qi cross section integrated over p"^ and made adimensional by 
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taking out a scale factor [1]. Correction to the naive parton model results comes en- 
tirely from the cross section which in QCD and in effective theories is quite different 
from the Q^-independent point-like cross section of the former. 

To analyze the scale dependence of the structure functions from (2.1), a conve- 
nient formalism is obtained by introducing the effective parton densities qi{x,Q^) 
and setting up evolution equations for them. It is well known that in perturbative 
QCD this is done in the leading logarithm approximation by resumming the ladder 
diagrams in the perturbation series, factoring it out and absorbing it in the defini- 
tion of the parton densities. In the case of the effective theories we are concerned 
with, it must be remembered however that (i) these theories are not asymptotically 
free, and (ii) the region of interest to which one would like to use these is not the 
asymptotic but the moderate region. 

All one can do here is expand the cross section in terms of the effective quark-pion 
coupling constant a^. 



a{z, Q') = d{z - 1) + ^a^{z, Q') + 0{al) (2) 

and redefine the parton densities as 

qi{x, Q'^) = qio{x) + Aqi{x, Q^) 

with 

Mi(x, Q') = ^J: [ %Ma'^(x/y, Q') + 0(al) (3) 
so that one recovers 

Ff{x,Q')/x = J2ehix,Q') (4) 

i 

which is the key formula in any theory of this kind. In (2.4), a'^J denotes the lowest 

order cross section for the process7*gj qiii. 

Then with the standard notation r = log{Q'^/ fi^) , //^ being a parameter defining 
the scale one can write down the evolution equation to the lowest order in 

^qi{x,T) =J2'Pij®Qj (5) 

where 

Vij ®qj= —Vij{x/y, T)qj{y, r) (6) 

y 

Vij{z,T)^^a'J{z,T) (7) 
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The last equation does not by itself guarantee universality or process-independence 
of the evolution kernels, an elegant feature of the splitting functions appearing in 
the leading log QCD evolution. 

The crucial point is that emission or absorption of an isovector causes a change 
of flavor of the initial quark. The lowest order cross section for the process 7*5 — > gV 
can be written as (after taking out the scale factor) 

ar (x, Q') = afix, Q') + ^^^(x, Q') (8) 

The cross-section for the process 7*7r — > q'q can be written as 

a^:^ = (x, Q') + - 1) Q') (9) 

where ad is the direct (t-channel exchange) cross section and cTj arises from the 
interference with the crossed channel (s-channel exchange) diagram. Both the terms 
are of the same order in and none of them can be neglected in the kinematic 
region of moderate and moderate energy to which the theory is intended to be 
applied. 

On the other hand, the cross section for the process j*q — > irq' can be written as 

ar(x,g^)=a7(x,g^) (10) 
This together with (2.8) leads to evolution kernels such that 

Vg',{l-X,Q^)j^V^gix,Q^) (11) 

This makes it impossible to sustain the probabilistic meaning of the evolution ker- 
nels. Further, the probabilistic interpretation would contradict the flavor- dependence 
of the evolution kernels Vun+ T^dir- as that would apparently violate the SU(2) 
symmetry. While this adds to the flavor-asymmetry of the sea distribution [9] which 
in any case would arise in such a theory due to a multiplicatively evolving part of the 
Gottfried sum [2] and may help explain the observed asymmetry, the probabilistic 
interpretation of the evolutions is clearly traded off. 



3 Evolution equation in chiral field theory 

For a speciflc case, we focus our attention to the chiral held theory [7] which was 
invoked by Eichten et al.[4] in their explanation of the NMC result on the Gottfried 
sum rule. The linearized effective lagrangian is 

Cu, = -^V^^^HtV^' (12) 



4 



Here 11 is the pseudoscalar octet of flavor SU(3) playing the role of Goldstone bosons, 
■i/' the quark field, qa the dimensionless axial vector coupling constant and / the pseu- 
doscalar decay constant. 



Contributions from direct cross sections 

The cross section for the process 7*? — >^ 7rg' arising from the direct (t channel) term 
is found to be of the form 

^d(7 q^nq) = — 2m, (m, + m,0 dt——-^-^ (13) 



where 



A' = min[-2ml + 2^{P^ + ml){K^ + ml) + 2PK, A^^^] (14) 
t2 = 2ml - + ^^)(^' + ml) + 2PK (15) 
P — m (16) 



^2 ^ js-ml- ml)'' - Amlml 
4s 

and is cut-off parameter of the chiral field theory and corresponds to the 

scale of chiral symmetry breaking. The fiavor factor has been suppressed in this and 
all the formulae given below. 

2. (Tdia*ii Qj^) 

The cross section for the process ~^ Q.j'^ arising from the direct (t channel) term 
is found to be of the from 

. * . (Jo 1 [^'^ ai + hit + cit' + dit^ 
aail - g.-) = j,^,^, dt ^^-^^ (18) 



where 



A" = mm[2^(p2 + rnl){K^ + ml) + 2PK - (m^ + m^), A^^s] (19) 
^2 = -2^(P2 + rn2)(ii:2 + ml) + 2PK + + (20) 
^2 ^ (s - - ml)'' - Amlml ^^^^ 



4s 

(g + g2_^2)2^4^2g2 

4s 



(22) 
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The expressions for the coefficients ai, 61, ci and di in the integrand are given in 
the appendix. 



3. 7*7r Qiqj 

The cross section for the process 7*7r q^qj arising from the direct (t channel) term 
is found to be of the from 

<T.(7 - - m,) -p^.J^. dt^^-^ (23) 



where 



= mtn[l{s + Q'-ml) + 2PK-mlAlsB] (24) 
1 

2' 



'2 



t2 = --{s + Q^ -ml) + 2PK + ml (25) 



4s 
s — Anig 



(27) 



The expressions for the ^2 and C2 in the integrand are given in the appendix. 
Contributions from interference cross sections 

1. ai{Yqi qjTi) 

The cross section for the process j*qi qjir arising from the cross term is found to 
be of the from. 

f * ^ 9a 1 1 f'^ a\ + b\t + (^ 

The expressions for the coefficients h\ and c\ in the integrand are given in the 
appendix. 

2. 7*7r qiqj 

The cross section for the process 7*7r qiqj arising from the crossed term is found 
to be of the form. 

/ * -N 9a ,^ 1 a'2 + bi^t + cit'^ 

<1 - - m,) = 7^32^(1 - -):^ "^' {t-mDit + al) ^^'^ 

where 

al = Qyz-{ml + ml) (30) 
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The expressions for the coefficients 03, 62 ^2 the integrand are given in the 
appendix. 

The cross section for 7*g — > nq' evaluated perturbatively in the lowest order 
leads to the evolution kernel V-n^q which has a z dependence 



dp^ 1 /"^^ , t 1 t2 dt2 

at- 



(31) 

where we have written = 9a/ P- Other quantities are defined as follows. 



= 2m' -J{P' + m2) {K^ + ml) + 2PK (32) 



ti = 2m5 - ^(p2 + ^2)(^2 + ^2) _ (33) 

where P and X are the magnitudes of the three momentum of the initial hadron 
and pion. In the above A^sb is an ultraviolet cut off parameter [4] which may be 
supposed to define the scale of chiral symmetry breaking. The evolution kernel 
Vq'q{z) corresponding to the process j*q — > q'l: is 

Vq,q{z)=gl{z,Q') + ^fnz,Q') (34) 
and the evolution kernel Vq^Triz) for the process 7*7r — >■ q'q is 

Vq>^{z) = g^Az, Q') + (if - 1) fi(^^ Q") (35) 

The expressions for g"^, g^, ff and are given in the Appendix. The quark densities 
are then found to satisfy the following equations 

duiz O'^) 

= Vuu®U + Vud®d+Vu^o®Tl° + Vun+ 7r+ (36) 



dr 

du{z, Q'^ 



dr 

and the pion densities satisfy 



Vuu®u + ®d + 'Pu-K- (8) 7r° + Vu-j:- ® (37) 
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^^^^^^:I^^V„+u®u + P,+a®d (38) 
dTr~{z, Q 



dr 

d7i°{z,Q^) 



= V^-d ^d + V^-u ® u (39) 

T.^nOg,^qi (40) 



dr 

The evolution kernels are explicitly flavor-dependent as stated in eqn(2.12). This 
implies that the Goldstone isotriplet has a non- vanishing contribution to the evo- 
lution of the non-singlet combinations q — u — d and q = u — d. This again is a 
consequence of inclusion of the interference term in the kinematic region of inter- 
est and is a feature characteristic of evolution through flavor-changing interactions. 
The moment equations resulting from these evolution equations will be discussed in 
section 5. The dependence of the evolution kernels on Q"^ is shown in Figures 1-5. 
The value of \sb is chosen to be llTOMeF. 
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Figure 1: Contribution of the t-channel diagram to Pq^riz) for 1. — 
2.0 GeV^soM line), 2. = 2.0 GeF^ (dashed hne), and 3. = 4.0 G'ey2(dotted 
hne) 

4 Probabilistic calculation of the evolution ker- 
nels 

In the kinematic region xP » Q > Pt, where P denotes the hadron momentum, 
the leading logarithmic approximation is known to be valid in QCD. The evolution 
kernels obtained in this approximation admits probabilistic interpretation and can in 
fact be calculated taking direct advantage of that. Durand and Putikka have shown 
that even the infrared divergences can be eliminated in such calculations if the 'loss 
terms' are correctly incorporated. In this section, we apply this method to derive 
the evolution equtions in the chiral field theory for the kinematic region specified 
above, the exercise leads to results which bring out striking differences between 
the evolution kernels in QCD and those in xFT. In particular, it illustrates the 
difficulties one would face if one tries to apply the latter to high energy. 

Following Durand and Putikka, we now denote by x the fraction of the hadron 
momentum carried by a parton. Using the vertex factor of x^T, the transition 
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Figure 2: Contribution of the interference term to the evolution kernel Pq 
1. g2 = 0.5 GeV^{so\id line), 2. ^ 2.O Ge^^ (dashed line), and 3. 
4.0 G'ey2(dotted hne). 
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Figure 4: Contribution of the interference term to the evolution kernel Pij{z). The 
plot does not show any perceptible variation with Q'^ in the range 0.5 — 4.0 GeV^. 
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probabilities can be calculated from these diagrams in a straightforward manner. 
Neglecting the quark masses and writing = 9a/ Pi the resulting evolution equa- 
tions are 



du 
d^ 



dTT+ 



(27r; 



+ 



j[4:xu{x) I dxi 

J o 

/ dxi[2xix^q~^{xi) + 

Jx 



xl{xi 



X] 



(xi — x)2 + e 



Xi 



X 



dr 



dr 



X 



1 dxi 
2x Jx xi 

1 dxi 



■UiXi, 



2(xi 



2xi 



X) 



-TT 



(+°) 



Xi 



X 



(xi — x)^ + e 



g+(xi) 

X [Xi 



{Xi 



xy + e 
- xY 



TT 



dxi — 
[x 



X 



— TT 



X] 



dxi — 
ix 



X 



Xx, 



(41) 



(42) 



(43) 



where e = p\/2P'^, q+ = u + d and 7r(+°) = 7r+ + 7r°. 

There are obvious differences between these equations and those obtained in 
QCD by the same procedure. The most important difference can be seen if one 
considers the limit e — > 0. Whereas inclusion of loss terms leads to elimination 
of infra-red divergences in the splitting functions in the case of QCD evolution, 
the divergence in the e ^ limit is not eliminated in the xFT evolution. This 
shows that the kinematic region corresponding to this limit is beyond the domain of 
applicability of the chiral field theory which is not valid above the chiral symmetry 
breaking scale. This is true for all such theories which would incorporate evolution 
through emission of chiral Goldstone bosons. 



5 Evolution of moments 

The Q^-evolution of the moments of the structure functions in the chiral field the- 
ory follows directly from the basic equations (36) - (40). However the resulting 
evolution is not as simple as in the case of the leading logarithmic approximation. 
Complication arises due to two distinctive features of the basic evolution equations. 
Firstly, the kernels of evolution arc Q^-dcpcndcnt. Secondly, the flavor-dependence 
of the pionic contribution makes it harder to decouple the moment equations. In 
particular, simple non-singlet combinations like q — q no longer evolve in purely 
multiplicative manner. 

We introduce the moments of the parton densities and the evolution kernel by 

Qn{T)= [' dx X^-'qr,{x,T) (44) 
^0 
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Tn^^'lr) = /' dxx^-^VAB{x, t) (45) 
Jo 

and take moments of the equations (36) - (40) to obtain 

^C/n(T) - 7n""t^n + 7n"'^n + 7n""°n„° + 7„"-+n„+ (46) 

|:f/n(T) = 7/^C/n + 7n^'5n + 7n^"°n°+7/-"n„- (47) 

^n„+(T) = 7n"^"C/n + 7n"^'^n (48) 

^n/(T) = i:7n'^°»'g^n (49) 

and similar equations for D„, Z)„, and n„~. 

Wc shall consider below the exact solution to these equations first and then make 
an approximation to bring out clearly the role of the pion decomposition process in 
the evolution of the moments. We shall omit the index n from the symbols Qn, 7^^ 
etc. in the expressions given below for reasons of typographical simplicity. 

5.1 Exact solutions 

Define the combinations U^~^ — U — U, n^~^ = n+ — 11" etc. and transform to 

- exp j^" dr7"") U^-\ 



Further writing 



X = 

a formal solution is obtained as 



X{t) = exp (^j^ M dr^ X(ro) (50) 



where 

M 



~du 



~ud 


^UTT 





~dn 


~TTd 






V r 

The evolution takes a simple form for the eigenfunctions of the matrix M. Therefore 
the first task is to diagonalize it. After diagonalization one finds that there are three 
eigenfunctions of evolution given by 

Vi{T) = + Tjd^-^ + Cn(-) (51) 
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where 



pTTU p 



^.pdu _j_ pdwpTTU 



A,; Af + T'^'^r'^'i 

The Aj's, being the three eigenvalues of M, are the roots of the cubic equation 

X^-aX + P = (52) 

Q, ^ ;j)«7r^^« _ y^^TTd _ ~du~ud ^53^ 
M — 77 7 — 'y'y')' (^54j 

The evolution of each such combination ^^(t) is then given by 

Viir) = Vi{Ti) exp(Ai(T, ti) 

It may be noted that the quantities F*-' are simply related to the moments of the 
cross section. Thus, putting back the moment index, 



5.2 The role of the pion splitting function 

Influnce of the pion splitting kernel on the evolution of the non-singlet first mo- 
ment can be clearly shown in a simplified approximate version of the moment 
equations(5.81)-(5.84). Let us introduce the symbols 7^^ = T^qiQj: T^q<^ — 'PqiQj 
j), K = Kq.^ Qm = K-K''^ Qm = K-K"" ^^^^re we assume the validity of 
the isospin SU(2) symmetry and charge conjugation invariance.The moments of the 
kernels defined above will be represented by 7^ etc, in obvious analogy with(5.80). 
Further, if we take the combination 

^{u-d) + {u- d) 
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and neglect the flavour-dependence of all other kernels except that of Vq-j^, the evo- 
lution equation for the combination of the quark and antiquark densities takes 
the simple form. 

-^g+(^,Q2)^g^^0g+ + p;^(7r+-7r-) (55) 
otr 

where V'^ — 'Pu-k+ ~ ^d7r+ ■ Now Qgg may be neglected as it involves a process which 
is second order in a^^. The equation for the first moment is then given by 

^Qt = K - 7,7)^^ + 7;ni (56) 

This equation can be solved to yield 

Qt{T) = A(t - To)[g+(To) + rdr'7;(r')ni(rOA(r',ro)] (57) 

J To 

where the evolution factor is given by 

A(r, To) = exp fb^g^r') - 7,7(^0]^^' (58) 

J To 

This also turns out to be 

A(t, To) = A expiU,^" (r) - ai^^ (r)] (59) 

A = expld^'^r^) - iai"°(ro)] (60) 

In the low region, the pion density may be supposed to be small enough for one 
to neglect the second term in (55) in comparison with the first term. This situation 
corresponds to a multiphcative evolution of the Gottfried sum which evolves from the 
quark model value to a lower value at moderate Q^, as observed experimentally [2] . 
At higher Q^, the evolution becomes more complicated due to the presence of the 
second term. 



6 Conclusion 

One can obtain Q^-evolution of the moments of the structure functions in the chiral 
field theory at low Q^, and the Gottfried sum does show variation with Q^. However, 
the kinematic region in which the theory is applicable makes it difficult to compare 
the results of the theory with the NMC and the other experimental results. Further, 
unlike the splitting functions in the Altarelli-Parisi formalisim, the evolution kernels 
in such effective theories admit simple probabilistic interpretation. 
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7 Appendix 



Given below are the coefficients appearing in the different cross sections in section 
3. Coefficients of the cross section for the process ^*qi qjir 



ai = mj(s + g") - 4m^m^g' - mj(4mj - 3m^) (61) 

bi = -2ml(s + Q^)+3ml(3ml + 2ml) (62) 

ci = s + Q^ — drUq — ml. (63) 

di = 1 (64) 

Coefficients of the cross section for the process 7*7r q^qj 

02 = mj(s + 2m5) +4mJm^(Q2 _^2^ (65) 

b2 = -2ml{s + 2ml) - Amlrnl (^^) 

C2 = s + 2ml (67) 

Coefficients of the cross section for the crossed channel process 'y*qi — > qjn 

^1 = ml{Q^ + 2ml + 3ml){s - nig - ml) - m^iml - ml){s + Sm^ - ml) 

+mg((5^ + "^g)(s + "^1 - m^) — 2ml{s{s — 2m^ - ml) + ml{m1 - ml)) 

-Am^mlis - 2m\ + ml) (68) 

6i = -2{s-m^^{s + Q'^ -m\-ml)+^m\{2m\ + ml) (69) 

4 = -2(s + m^) (70) 

Coefficients of the cross section for the crossed channel process 7*7r — > qiqj 
4 = -mj(g2 + m2)(s-2mj)-m;(g2 + ^2)(s + 2g2_2m2)-2m2(s + g2-m: 



2\2 



+2mlml{s + g^ + m^) - 4m^m2 (2m^ + Sm^) + 8m^m^(s + g^) (71) 
6*2 = -{s -2ml){2s + Q^ -2ml-2ml) ~ s{Q^ + ml) -ml{s + 2Q^ -2ml) 

-2m\{2s + 2g2 - 2ml - + ^rn^ml (72) 

4 = -2s (73) 
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